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ABSTRACT 

This thesis describes the results of a theoretical study 
to predict the thermal behavior of an internally-cooled tube 
in a condensing vapor. 

The analysis constitutes a unique application of the 
finite element method and provides new insights into the 
effects of circumferential conduction upon condenser tube 
performance. Comparisons are made between the present analysis 
and the theoretical and experimental works of others. The 
inclusion of circumferential conduction leads to an improve- 
ment in the predictive capabilities of the analytical model. 
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I. INTRODUCTION 



A. BACKGROUND INFORMATION 

As the size, complexity, and hence the cost of modern 
warships have grown in the last thirty years, advances in 
basic research have been made and the computer has become a 
useful tool to the researcher and designer. The application 
of up-to-date research and computer technology to the design 
of modern warships is common but there remain major, costly 
components relatively unchanged and unreviewed in decades. 

It is in the interest of the Navy to ferret out these compo- 
nents and re-examine them. 

One such component is the naval condenser. The basic 
engineering developments currently used to design naval con- 
densers are from the Heat Exchanger Institute (H.E.I.) 
Standards for Steam Surface Condensers [1] and the standards 
of the Tubular Exchanger Manufacturers Association (T.E.M.A.) 
[2]. The Navy has further utilized the H.E.I. information 
in the Design Data Sheet (D.D.S.) [3]. While these standards 
have proven reliable, their performance margins and extensive 
use of averaged values may result in a larger, costlier con- 
denser at no increase in performance and with an unnecessary 
margin of reliability. These performance margins, once 
necessary to insure reliability, are now called into question 
as the mysteries of condensation are unravelled. 
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It is for the purpose of reducing the size of the con- 
denser that many researchers are trying to increase the heat 
transfer rate per unit of condensing surface area. 

B. CONDENSATION IN A NAVAL CONDENSER 

A typical naval condenser uses the heat transfer mecha- 
nism of filmwise condensation to the outside of horizontal, 
cylindrical tubes. Dropwise condensation is therefore not 
considered in this thesis. 

The inside of the tube has single-phase turbulent forced 
convection heat transfer. The correlations associated with 
this phenomenon are accurate enough to justify confidence in 
their use. The fouling of the inside of the tube is poorly 
understood, and varies greatly, so that a necessary design 
allowance is made and a performance margin created. The 
heat transfer through the tube wall itself is by conduction, 
a classical problem in cylindrical coordinates, and can be 
calculated as accurately as the properties of the tube material 
and the surrounding heat transfer coefficients are known. 

Aside from a greater knowledge of coolant-side fouling 
and how to control it, it is largely on the steam side of the 
tube that advances of knowledge can be reasonably expected. 

Here there are three important phenomena. The first is the 
transfer of heat through a thin layer of condensate from 
a vapor, either quiescent or flowing, which causes that vapor 
to lose its latent heat of vaporization and join the film of 
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condensate on the tube. As condensing particles join this 
liquid layer they retain some, none, or all of the momentum 
they possessed as a vapor. The dynamic interaction of the 
vapor and the liquid under both forced flow and quiescent 
vapor conditions presents a formidable problem to the 
researcher. A second phenomenon on the steam side of the 
tube is the presence of noncondensible gases which complicates 
and usually degrades the performance of the condenser. The 
third phenomenon is the inundation of a condenser tube by 
"condensate rain," or condensate condensed elsewhere in the 
condenser which, as it makes its way to the hotwell, impinges 
upon the condenser tube. 

C. IMPROVING THE HEAT TRANSFER 

If the rate of heat transfer for a given area of condensing 
surface were to be increased, then the total size of a con- 
denser could be reduced. The result would be savings in space, 
weight, and possibly cost. Various methods of increasing 
the heat transfer in a naval condenser have been proposed. 

Among these are: 

1 . Reducing the Internal Fouling 

The reduction of internal fouling would certainly 
increase the rate of heat transfer. While various tube metals 
and chemical treatments to the coolant side of the condenser 
can reduce fouling, operational and other considerations 
have prevented their use in naval condensers. Aside from 
mechanical cleaning, no acceptable method of significantly 
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reducing marine fouling in a naval condenser has presented 
i ts el f . 

2 . Eliminating Noncondensible Gases 

The elimination or reduction of noncondensible gases 
in the condensing steam would be expected to increase the 
heat transfer rate. It is a fact of realistic operation in 
a marine environment that noncondensible gases cannot be 
eliminated. There is room for improvement in their continual 
removal from the condenser by air ejector devices. As long 
as complicated systems with mechanical seals are used in the 
shipboard environment, some amount of noncondensible gases 
will be in the condenser. 

3 . Enhancing Heat Transfer 

a. Shaping of Tubes 

Various ingenious shapes of tubes have been de- 
monstrated to increase heat transfer rate for a given amount 
of tube surface area. Complications arise from these tubes, 
however, not the least being a sharp increase in the pumping 
power required to maintain an adequate coolant flow rate. 

b. Promoting Dropwise Condensation 

The use of various promoters upon the tube outer 
surface to prevent the forming of a liquid film, thus keeping 
the condensate in droplet form, has been proposed and demon- 
strated by many authors. Since it is not the intention of 
this thesis to examine the effect of this phenomenon, suffice 
it to say that until now there remains a question as to whether 
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any promoter has demonstrated the robustness to withstand 
continued, long-term use. 

4 . Eliminating or Reducing Inundation 

Since condensers are formed with banks of tubes, it 
is inevitable that condensate draining from one tube should 
impinge upon some other lower tube. The effect of this in 
filmwise condensation is normally to thicken the condensate 
film and to reduce the rate of heat transfer. The proper 
placement of collection trays or baffles is of great impor- 
tance in channelling falling condensate away from lower tubes. 
The number and placement of these collection trays is not 
examined in this thesis. 

5 . Utilizing High Vapor Velocities 

In many regions of a naval condenser, vapor velocities 
of 300 ft/sec or more may occur. These high velocities create 
a forced flow situation on the steam side, thin the film of 
condensate, and generally increase the heat transfer coeffi- 
cients on the steam side. The higher velocities also cause 
complications, one of which is the separation of the vapor 
boundary layer causing the condensate film to thicken suddenly. 
The heat transfer in the region of the tube which is after 
the vapor phase separation point is relatively unknown. 

D. REVIEW OF GENERAL RESEARCH IN CONDENSATION 

In analyzing the condensation of steam on tubes, the 
cornerstone of research was laid by Nusselt [4] who analyzed 
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a quiescent vapor condensing on a vertical flat plate and on 
a horizontal cylinder. Nusselt's analysis allowed for no 
interfacial shear stress between the flowing condensate and 
the condensing vapor. He also used a constant tube wall sur- 
face temperature. 

Rohsenow, et al . , [5] later expanded this work to include 
the effect of vapor shear at the vapor-liquid interface upon 
the condensation rates. This analysis was for a vertical 
flat plate of uniform temperature with either turbulent or 
laminar flow of condensate. 

Sparrow and Gregg [6] conducted the first analysis of sur- 
face condensation using boundary layer theory. They included 
momentum terms in the condensate film analysis to account for 
acceleration of the flowing condensate. They showed that the 
momentum terms were insignificant under most conditions. 

Koh, et al . , [7] did further work allowing the falling 
condensate to drag the vapor along with it, while still using 
a constant wall surface temperature. 

Shekriladze and Gomelauri [8] expanded upon the earlier 
work of Rohsenow, et al., [5], analyzing the effect of a 
flowing vapor upon laminar filmwise condensation in the 
absence of significant gravity effects. Their model used a 
constant wall surface temperature. 

Denney and Mills [9] expanded this area of analysis to 
include the influence of gravity upon the condensate film, 
holding to a constant wall surface temperature. 
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E. RECENT WORK CONCERNING CONDENSATION ON TUBES 
1 . Vapor-Condensate Interfacial Shear Models 

One of the areas of continual interest in condensation 
of steam on a tube is the question of the effect of the shear 
or drag of the condensing vapor upon the flowing condensate. 
Theories which attempt to quantify the interaction abound, but 
recent ones take one of two general forms. 

The first type of model says that total shear is the 
result of the summation of dry friction shear of the vapor 
upon a bare tube and a shear caused by the condensing vapor 
retaining some, none, or all of its vapor phase momentum. The 
deceleration of this vapor as it joins the liquid layer causes 
this shear. Various authors [8, 10, 11, 12, 13,] subscribe 
to some combination of these two shears. 

Other researchers [14, 15] subscribe to a theory which 
reaches more deeply into the molecular level interactions and 
in reality appears to be a more minute examination of the 
same idea. Using a Reynolds Flux concept, these authors state 
that for a flowing fluid there is a transport of momentum in 
a direction normal to the mass flow which is on the molecular 
level in laminar flow but involves larger amounts of fluid in 
turbulent flow. This transport is a function of temperature 
in laminar flow and of temperature and free stream turbulence 
in turbulent flow. As one of these cross-flowing particles 
approaches a solid wall, it is decelerated to a standstill. 

As it rebounds from the wall, it is reaccelerated. Thus, 



an overall rate of condensation can be represented by super- 
inposing a velocity perpendicular to the wall, similar to 
suction through a dry wall. 

The latter theory includes within it a relationship 
between the rate of condensation and temperature drop that 
must necessarily occur in the process [14]. The former 
frictional shear force theories commonly assume that the 
outer surface of the condensate film is at the steam 
saturation temperature. 

2 . Wall Temperature Models 

Fujii, et al . , [16] display from their experimental 
data that tube temperature varies around the outside surface 
of the tube. In their calculations, however, they use averaged 
values for tube wall temperature and for the temperature-depen- 
dent heat transfer coefficients to calculate heat transfer.^ 
Nobbs and Mayhew [17] show that for steam flowing 

vertically downward, a variable temperature analysis results 
in a marked overall change in heat transfer. The state that 
earlier isothermal tube wall surface models were optimistic 
in their predicted heat transfer. 



Since the writing of this theis, new results have been 
obtained by Fujii and his coworkers . The reader is referred 
to the Proceedings of the Workshop on Modern Developments in 
Mari ne Condensers , Naval Postgraduate School, 26-28 March, 1 980. 



Nicol and Wallace [13] show that local heat transfer 
coefficients are indeed a function of local tube surface 
temperature and recognize the importance of the circumferential 
path in heat conduction. In one of the first analyses to 
include this path of heat transfer, these authors apply a 
relaxation technique and iterative scheme, changing the 
temperature distribution and then the local heat transfer 
coefficients. It is interesting to note that their results 
show that when the circumferential heat conduction path is 
included, the average heat transfer coefficient rises about 
twenty-five percent. 

Nicol, et al . , [18] analyze steam crossflow upon condenser 
tube (Fig. 1) using a formulation wherein local tube surface 
temperatures are used to calculate local heat transfer co- 
efficients. Thier results show that the variable wall 
temperature formulation predicts less heat transfer then older 
isothermal wall models. The effect of circumferential con- 
duction around the tube wall is not included in their analysis. 
A close look at Fig. 7 of their work, partially reproduced here 
as Fig. 2, shows a theoretical wall temperature profile and 
a few measurements of surface temperature around the tube wall. 
In their discussion of this particular figure they state: 

"The profiles are seen to agree quite well, the main 
di ferrences being in the rounding off of the more extreme 
peaks of the theoretical profile because of the circumferential 
conduction in the tube wall, which is of course more pronounced 
in those regions where the temperature gradient is high." 
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Nicol, et al . , [18] also present theoretical curves of 
local heat transfer coefficients, local heat flux, and local 
tube surface temperature. These cases are presented for both 
an older isothermal wall model and an anisothermal wall model 
where heat can flow only radially in the tube wall. 

Fujii, et al., [19] use the results of Nicol, et al., 
[18] for variable wall surface temperature, coupled with 
earlier work [20] and develop another model based upon the 
assumption of uniform heat flux around the tube. These 
authors also recognize the path of ci rcumferenti al heat con- 
duction as an important area for further work. 

F . Purpose of Research 

It is in this final area that the author has analyzed, 
with the computer, a single tube of a condenser and has 
included in his analysis the circumferential conduction of 
heat in the tube wall, a phenomenon not customarily included 
in condenser heat transfer analysis. 
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II. A SPECIFIC PROBLEM FOR STUDY 



A. THE SINGLE TUBE PROBLEM 

A review of prior research has indicated ideas and trends. 
One basic idea was that local heat transfer coefficients de- 
pend upon local temperature differences and that the shearing 
force depends upon these heat transfer coefficients. The 
trend was away from older isothermal wall models to newer 
variable wall temperature ( ani s othermal ) models. The logical 
continuation of the trend would indicate that a wall fully 
capable of two dimensional heat conduction would have been 
be t ter still. 

Since actual tube walls can conduct heat in the circum- 
ferential direction, then a model which could account for 
this would be an improvement over previous models. In order 
to test the plethora of shear force and condensation rate 
theories, a model including the effect of circumferential 
conduction is desirable, and may prove to be essential. 

B. THE PROBLEM STATEMENT 

The problem was posed for the upper half of a cylindrical 
condenser tube with a cross-flowing vapor, that is, where the 
direction of the gravitational force and the direction of 
vapor flow where perpendicular, (See Fig. la). Symmetry with 
the other half of the tube was assumed, allowing the simplifi- 
cation that the ends of the half tube, where it would normally 



join the other half, were adiabatic walls. The problem for 
solution was set up so that e, the local temperature differ- 
ence, was defined thusly: 



0 (r,«j> ) = T - T ( r , 4 > ) 



( 1 ) 



While T s may vary with pressure changes around the tube, it 
is assumed constant herein. The governing equation for heat 
transfer within the tube wall is the Laplace Equation: 
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with the following boundary conditions: 

(fr) r. .* = h i [(J<r i>* ) - e c ] 
= h.e. 



= k w 



% = - k w 



m 



r„ ,+ 



I 0 



On the inside of the tube turbulent forced convection was 
assumed and the Di ttus-Boel ter correlation [21] was used to 
provide a heat transfer coefficient. 

On the outside of the tube, the condensate flow was taken 
to be laminar, and heat transfer through the thin condensate 
layer was assumed to be purely by conduction. For such con- 
ditions the film heat transfer coefficient is given by: 



\ -r 



(3) 



The film thickness 5 was calculated from the formulation due 
to Nicol, et al . , [18] (see also Appendix A). This analysis 
yields the following first-order differential equation: 
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The shear term, t , used by Nicol, et al., [18] was one 
given by Schlicting [22] and is of the following form: 
p U “ 

V V ,( 6.973$ - 2 . 732$ 3 + 0 . 29 2$ 5 - 0.01 83$ 7 
u v r o p v \ ^ 

p v ) + 0 . 000043$ 9 - 0 . 0001 1 5$ 1 1 ) . (5) 

For the sake of comparison this author used the same shear 
stress model. This model calculates vapor boundary layer 
separation at about 108 degrees from the forward stagnation 
point. Thereafter, the shear stress was set to zero. 

In order to integrate Eq. (4) a starting value of con- 
densate film thickness at the stagnation point was needed. 
Because high velocity steam flowing horizontally normal to a 

H £ 

horizontal tube was to be studied, was set equal to zero 
for the stagnation point {$ =0), and an approximate starting 
value for the condensate film thickness, 6, was found from 
Eqs . (4) and ( 5) . 

Once this starting value was in hand, Eq. (4) was inte- 
grated, in its entirety, around the tube using an estimate of 
the tube surface temperature difference, [e Q = 9 ( r Q , $ ) ] . The 
integration resulted in the local film thickness around the 
tube wall. Local heat transfer coefficients for the steam 
side were then calculated according to Eq. (3). 

Using these local heat transfer coefficients Eq. (2) could 
then be integrated, resulting in a temperature distribution 
throughout the solid tube wall. With a new set of tube-wall 
outer surface temperatures in hand, a new set of local heat 




transfer coefficients could be calculated and applied for 
another iteration of the Laplace equation. If temperatures 
converged and heat transfer balanced, then a solution was in 
hand . 

C. THE GALERKIN FINITE ELEMENT METHOD 

The purpose of this application of the Galerkin Finite 
Element method is to solve a steady state two dimensioanl 
heat conduction problem by converting the governing partial 
differential equation, the Laplace equations in this case, 
to a series of algebraic equations. Each of these algebraic 
equations represents the partial differential equation over 
a samll sub-region of the total domain called an element. 

The result of the simultaneous solution of these algebraic 
equations is to drive the integral square error of the 
approximate solution to zero over each element. Elements may 
be of different sizes and may be clustered in areas where 
greater accuracy is desired. 

The Galerkin Finite Element method differs from the 
Galerkin method in that, since the algebraic equation repre- 
senting the approximate solution exists over only one small 
element, the resultant matrix of algebraic equations for 
solution takes on a banded form. That is, it has zero upper 
and lower triangles of significant size. The resultant 
savings in computer storage space and computational time by 
using special matrix equation solvers is significant. 



In this discreet approximation method the boundary con- 
ditions are applied locally to each element on the boundary 
of the domain. In this case linear triangular elements 
were used. The algebraic approximations to the solution, 
called basis or shape functions, were linear. The resultant 
approximate solution amounts to approximating a smooth three 
dimensional surface by a mosaic containing a large number of 
small triangles. Triangles size varied, small ones being 
used in such areas of interest as the forward stagnation 
point and anywhere a steep temperature gradient could be 
expected to exist. Larger elements were used in areas where 
the temperature profile was considered to be relatively 
unchanging. As solutions progressed, refinement of the 
triangular meshes highlighted and concentrated upon areas of 
i nteres t . 

For greater detail on both the Galerkin method and the 
Galerkin finite element method, the interested reader is 
referred to Ozisik [23] and Fairweather [24]. 

The formulation and computer program were tested using 
classical two-dimensional problems of heat conduction from 
Carslaw and Jaeger [25]. The first of these was a rectangle 
with one side at a given temperature and the other three at 
a specified lower temperature. The second was a rectangle 
with two adjoing sides adiabatic, a third side at a given 
temperature, and the fourth having a given heat transfer 



coefficient to a medium at some lower temperature. In both 
cases the finite element solutions converged to exact solutions. 

A further test of the program and formulation was the 
problem of a rectangle with two opposite sides adiabatic, a 
theird side at a given temperature, and the fourth side having 
a local heat transfer coefficient to some medium at a lower 
temperature. This local heat transfer coefficient was 
dependent upon both the local surface temperature and the 
position along the wall. Since this problem cannot be compared 
to a classical solution, any solution could only be given 
credence by detailed examination. It was seen that the 
solution converged for smaller and smaller triangle sizes, 
that plots of isotherms internal to the rectangle were smooth 
and not obviously in error, and that heat transfer into and 
out of the rectangle was in balance. While certainly no 
proof, these comparisons led to confidence in the program. 

D. APPLICATION OF THE INNER BOUNDARY CONDITION 

For the inside of the tube, the boundary condition of con- 
vection heat transfer, Eq. (2), can be simple represented by 
noting that the inside film coefficient, h. , is given by the 
Dittus-Boelter correlation [21], and is a function of coolant 
bulk temperature and flow rate. Thus, since the thermal 
resistance at his interface is independent of local tube wall 
temperature, the resistance may be represented by a fictitious 
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thickness of insulating mate rial according to the relationship: 



R = 



1 



r i 1n iB - tr^ 



h i A 



K f A 



ri - tf 



( 6 ) 



From this, the thickness of the fictitious ring t^, was deter- 
mined and the equivalent boundary condition, that the coolant 
temperature existed at the inside of the fictitious ring, was 
i mposed . 



E. SOLUTION PROCEDURE 

On the outside of the tube, the local heat transfer co- 
efficients were computed all around the tube for the current 
value of the calculated local outside wall temperature. These 
local heat transfer coefficients were applied to elements 
having sides on the outer tube wall. 

In order to converge to the proper wall temperatures and 
temperature-dependent heat transfer coefficients, an iterative 
scheme was followed wherein a new set of local heat transfer 
coefficients for the outside of the tube were computed from 
new solution temperatures. The process was repeated until 
temperatures at all nodal points changed less than 0.02 degrees 
Celsius for consecutive iterations. 

Having stabilized the local temperatures throughout the 
tube wall, a check was made to insure that the heat transfer 
balanced. The total heat deposited upon the condenser tube 
outer wall by the condensing steam was compared to the total 
heat conducted into the tube wall at the outer surface and to 
the total heat conducted from the tube at its inner surface. 
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F. COMPARISON WITH PREVIOUS SOLUTIONS 

The aforementioned paper of Nicol, et al . , [18] provided 
fertile ground for comparison, but only if the test conditions 
could be simulated. These authors were able to achieve a 
given average tube outer surface temperature by varying the 
flowrate of the coolant and hence the internal heat transfer 
coefficient. In the finite element formulation this condi- 
tion is obtained by varying the thickness of the aforementioned 
fictitious ring. In order to preclude significant circum- 
ferential conduction in the fictitious ring, the circumferen- 
tial thermal conductivity of the ficticious ring was chosen 
to be one one-hundredth of the thermal conductivity of the 
tube wall. 
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III. RESULTS 



A. TUBE WALL TEMPERATURE 

The tube wall temperature predictions of the author's 
finite element model are shown in Fig. (3), and are seen to 
support the prediction of Nicol et al . , [18] that circum- 
ferential heat conduction would substantially round off the 
theoretical temperature profile obtained by considering 
radial heat conduction only. It is also worthwhile to note 
that the maximum and minimum temperatures are less extreme 
than those predicted by the radial heat flow model of Nicol 
et al., [18] and generally agree more closely with their 
experimental results. While the actual differences in 
temperature appear small, near the separation point of the 
vapor flow these differences approach thirty percent. 

Nicol et al . , [18] also show a theoretical tube wall 
surface temperature profile for another average wall tempera- 
ture. This profile is reproduced in Fig. 4 along with a 
theoretical profile developed from the author's model. As 
before, the steep gradients are rounded off and the overall 
temperature profile is less extreme. 

A plot of the isotherms within the tube wall is quite 
revealing and is shown in Fig. (5). The direction of vapor 
flow and point of vapor boundary layer separation, as indicated 
by the thickening of the condensate layer, are shown on the 
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figure. The isotherms are seen to start close together near 
the vapor stagnation point and to spread, reflecting lower 
rates of heat transfer, progressively around the tube. In 
the vicinity of separation of the vapor boundary layer, the 
isotherms are seen to deviate significantly from their con- 
centric upstream orientation. This reflects a significant 
amount of heat transfer circumferentially around the tube in 
this general vicinity. The much wider spacing of the isotherms 
on the back part of the tube indicates that according to this 
model it is an area of relatively low heat transfer rate. 

1 . Local Heat Transfer Coefficients 

Figure 6 shows a comparison of local heat transfer 
coefficients around the tube wall. The isothermal wall and 
radial heat flow local transfer coefficients of Ni col et al . , 
are partially reproduced. The results of the finite element 
model are seen to lie generally between the other two models 
with the exception of the vapor-separated region. In the 
region prior to separation, as seen in Fig. 4, the wall 
temperature difference of the finite element model lies be- 
tween the isothermal wall case and radial heat flow case. 

From Eq . (4) at the stagnation point, the initial value of 
the film thickness is proportional to 9 Q 1 / 3 . Therefore, a 
larger temperature difference will create a larger initial 
film thickness, and therefore a smaller heat transfer 
coefficient, as given by Eq . (3). It is not unexpected that 
this formulation and solution should predict a greater heat 
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transfer coefficient in the separated region than the iso- 
thermal wall model. Since there was a higher wall temperature 
and therefore less steam condensed in the unseparated portion 
of the tube in the present model, there is a thinner conden- 
sate layer in the after section of the tube. In the separated 
region the greater heat transfer coefficient predicted by 
the author's model over the radial heat flow model is not 
fully explained, but is no-doubt due, in part, to the assumption 
of symmetry in the present formulation. 

2 . Local Heat Flux 

In still a third examination under these same condi- 
tions, Nicol et al . in their Fig. 6 present local heat flux 
for both the isothermal wall and radial heat flow models. 

These profiles are partially reproduced here in Fig. 7 along 
with the author's results which again, as expected, lie between 
the two previous models except in the separated- vapor region 
of the tube. It is interesting to note that the inclusion 
of circumferential effects generates a heat flux which is 
much closer to being uniform than that which is obtained with 
the isothermal wall model. 
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IV. CONCLUSIONS 



The circumferential path of heat conduction has been 
successfully included in the solution of laminar filmwise 
condensation of high velocity steam on condenser tubes. The 
author's anisothermal fully-conducting wall model yields 
closer agreement to previously measured tube wall surface 
temperatures than a previous anisothermal wall model which 
allowed only radial heat flow. The maximum and minimum values 
of wall temperature were less extreme for the anisothermal 
fully conducting wall than for the anisothermal wall with 
only radial heat flow. 

The addition of a circumferential path for heat flow 
results in a complicated interplay between this heat flow 
and heat being deposited locally by condensation and leaving 
in the internal cooling water. The addition of vapor velocity 
effects on the outside surface of the film, and especially 
the sudden changes arising from separation of the vapor flow 
result in a highly complex sequence of thermal and fluid- 
dynamical interactions. The methods developed herein, 
though s ti 1 1 s ub ject to comprehensive experimental verification, 
are adequate to predict these complex phenomena. 
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V. RECOMMENDATIONS 



Now that a model for two-dimensional heat conduction which 
includes temperature-dependent local heat transfer coefficients 
has been developed and demonstrated, the use of this model for 
further study is suggested. Areas thought fruitful for further 
examination are both the shear stress models and the tempera- 
ture-dependent condensation rate models. For shear stress 
models the separated vapor region or trailing area must be 
given closer examination since it accounts for about forty 
percent of the surface area of the tube and the fluid flow 
characteristics in this region remain relatively unknown. In 
the study of temperature-dependent condensation rate models, 
further examination of the effect of the velocity of the 
flowing vapor and local temperature is expected to be fruitful. 

The vigorous exercise of this model through a broad range 
of steam-side and coolant side conditions is suggested. 

The model of condensation with constant heat flux pro- 
posed by several authors [18, 20, 26], should be re-examined 
since circumferential heat conduction in the tube wall, when 
included in theanalysis, seems to yield a nearly constant 
heat flux distribution. 
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APPENDIX A 



Local Rate of Change of Condensate Film Thickness 

A force balance on the shaded element of Fig. lb yields 
the following: 

T v r o d(|> " T £ r o d * ' 9 (p a -p v )cos<|> r Q (<5-y) d<j> = 0 



so that: 



3 U 



= t - (p.-pjg c o s <j> (<S-y) 



sl = u 3y V v H v 



Integrati ng : 

u = _.Xy _ (pjj_-p v )^- cos4> ( <S y-^- ) 
The average velocity, u is: 

; = 1 / ,, d, = i L f . 2 ( 



6 "“o 



6 u 



U 



COSij) 



T V^ /p p \q <$ 2 

2T ' ( * v) 7 cos * 3 



i 3 i 3 

2 ' 6 



, or 

(8) 



The mass flow rate of condensate around the tube is: 



rii„ = p .5 u 

and , wi th Eq . ( 8) , 



„ T v P 
\ -% 



V V 



gp t( p r°ii) 

3y 



cos<J> 6 3 



The rate of change of mass flow rate with respect to angle is 



dl \ _ p Jl' S dT v + p 2 t v <S d<$ + gp Ji( p r p v)5 3 s 1 n4> 
d<J> 2y d<J> y d(j> 3 y 

9P £( P g,- P v) COStf 6 2 dS_ 
u d(f> 
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Since drf)/d<f> can also be equated to the amount of steam 



condensing, the following relationship applies 



drfij 
d <f> 



r o k * 9 o 



<$ h* 



fg 



where h 1 ^ is the latent heat of vaporization corrected for 
subcooling. By making this substitution and solving for 
d <S / d <j> , the relationship given in Eq. (4) is obtained. 
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APPENDIX B: PROGRAM LISTING 
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APPENDIX C 





• (a) 




(b) 



Figure 1 . Schematic diagrams of physical model 
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Outside wall temperature, 




Figure 2. Predicted and measured T q versus angle from stagnation point [18]. 
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Figure 3. Predicted and measured T versus angle from forward stagnation 
point. 
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Figure 4. Predicted T Q versus angle from forward stagnation point. 
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Figure 5. Isotherms (In degrees Celsius) of 6 Inside tube wall for top 
half of tube. Thickness of condensate film Is shown with a 
magnification of 25 times. 



Local heat transfer coefficient, h», (W/m 2 o C)xlO 




Figure 6. Predicted local heat transfer coefficient versus angle from 
forward stagnation point. 
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Figure 7. Predicted local heat flux versus angle from forward stagnation 
point. 
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